Induced scalarization in boson stars and scalar gravitational radiation 
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The dynamical evolution of boson stars in scalar-tensor theories of gravity is considered in the 
physical (Jordan) frame. We focus on the study of spontaneous and induced scalarization, for which 
we take as initial data configurations on the well-known S-branch of a single boson star in general 
relativity. We show that during the scalarization process a strong emission of scalar radiation occurs. 
The new stable configurations (S-branch) of a single boson star within a particular scalar-tensor 
theory are also presented. 
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I. INTRODUCTION 

Scalar-tensor theories (STT) are alternative theories of 
gravitation where a spin— degree of freedom <p can ac- 
company the usual tensor spin— 2 modes (see Ref. [l| for 
a review). There are two mathematical representations 
of the STT: 1) The physical frame (also known as the 
"Jordan" frame), where test particles follow geodesies of 
spacetime and the scalar field <f> couples non-minimally 
to the curvature: 2) The Einstein frame, obtained by 
a conformal transformation of the metric, where the 
scalar field couples minimally to the curvature and non- 
minimally to the matter fields 0. 

STT are perhaps the simplest, well motivated and 
most competitive theories of gravitation after General 
Relativity (GR), the most prominent example being the 
well-known Jordan-Brans-Dicke theory 0, [if . Intuitively, 
STT can be seen as theories with a varying effective gravi- 
tational "constant" . Although so far there is not observa- 
tional evidence that such scalar gravitational field exists, 
one can use the experimental and observational tests of 
GR toput limits on its existence and its possible interac- 
tions @. Using data from the binary pulsar, for instance, 
it is possible to put limits on some classes of STT which 
restrict the form of the non-minimal coupling (NMC) to 
the curvature. Nevertheless, these bounds still allow a 
NMC constant of order unity . 

Despite the fact that STT were proposed several 
decades ago, it has only been recently that several phe- 
nomena associated with them, and with no counterpart 
in GR, have been analyzed. For instance, in the cosmo- 
logical context, STT have been proposed as alternatives 
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to the cosmological constant in order to explain the ac- 
celerated expansion of the Universe [71412} . 

In the astrophysical scenario, Damour and Esposito- 
Farese @, [l~3| discovered that neutron star models within 
STT may undergo a phase transition that consists on the 
appearance of a non-trivial configuration of the scalar 
field </> in the absence of sources and with vanishing 
asymptotic value. This phenomenon has been named 
spontaneous scalarization (SS) due to its similarities with 
the spontaneous magnetization of ferromagnets at low 
temperatures. The stability analysis for the transition 
to SS was first performed by Harada [HI, EH] ■ It is now 
understood that SS arises under certain conditions where 
the appearance of a non-trivial scalar field gives rise to 
a stationary configuration that minimizes the energy of 
the star with fixed baryon number. 

Further analysis have confirmed that the SS phe- 
nomenon takes place in neutron stars independently of 
the equation of state (EOS) used to describe the nuclear 
matter (l6l-[T8j . In boson stars this phenomenon was first 
studied by Whinnett [l9|], who constructed stationary 
scalarized configurations with a self-interaction potential 
for the scalar field. More recently, the dynamic transition 
to SS was analyzed by us in the Jordan frame and with- 
out self- interaction [20| . One important feature of this 
phenomenon is that it can occur even when the parame- 
ters of the theory satisfy the stringent bounds imposed by 
the Solar System experiments, notably, when the Brans- 
Dicke parameter is chosen to be arbitrarily large. 

The SS phenomenon is accompanied by the "sud- 
den" appearance of a new global quantity termed scalar 
charge, where by sudden we mean that the derivative 
of this charge with respect to the central energy den- 
sity at the critical point is infinite. The scalar charge is 
the analogous of the magnetization of ferromagnets men- 
tioned above. Moreover, just as in the neutron star case, 
in boson stars one also finds that beyond a certain crit- 
ical central energy-density, the stationary configurations 
that are energetically preferred are those where the SS 
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ensues. 

A phenomenon similar to SS, but that occurs when 
a background scalar field is present, is called induced 
scalarization (IS). It corresponds to the case where the 
scalar field does not vanish asymptotically. In this situ- 
ation the scalar charge does not exhibit a discontinuous 
"jump" as the object becomes more compact, and the 
transition to scalarization is smoothed out by the pres- 
ence of the background field. 

Another important feature of STT is the prediction 
of scalar gravitational waves. While GR predicts only 
quadrupolc gravitational radiation in the "far zone" , 
STT predicts the existence of monopolar gravitational 
waves that can be emitted even in the case of spherical 
symmetry plj . The new polarization of this scalar mode 
is of breathing type since it affects all directions isotropi- 
cally 13. 

A simple scenario where such scalar gravitational 
waves might be produced is precisely during the scalar- 
ization process of a spherical compact object. The am- 
plitude of such waves (in the linear approximation) is 
linked directly with both the form of the NMC and the 
asymptotic value of the scalar field. For instance, when 
this asymptotic value vanishes, it turns out that certain 
classes of STT do not lead to the emission of scalar grav- 
itational waves. This implies in particular that in such 
classes of STT the SS phenomenon does not produce 
monopolar waves, rather it is only in the IS scenario that 
such theories can lead to an emission of scalar gravita- 
tional radiation. Since we will be working with one such 
class of STT, it is then particularly important to make a 
clear distinction between the SS and the IS phenomena. 

The first dynamical analysis of the scalarization phe- 
nomenon in neutron stars was made by Novak [jq . He 
did not only confirm the dynamical transition to the 
scalarization state, but also the emission of scalar grav- 
itational waves. Moreover, he also studied the emission 
of scalar gravitational radiation when a scalarized star 
collapses into a black hole [22j ■ 

Recent studies of neutron star oscillations within STT 
have shown that, in addition to the emission of scalar 
gravitational waves, the quadrupole gravitational radia- 
tion is also disturbed as compared to the corresponding 
signals in GR [23[ . Therefore, even if the detection chan- 
nels for scalar gravitational waves are "switched off" , the 
detection of gravitational waves of spin-2 coming from 
these sources might still validate STT or put even more 
stringent bounds on their parameters. Of course, the 
direct detection of scalar gravitational waves (or the ab- 
sence thereof) would also help to discriminate between 
several alternative theories. 

In this work, we present a systematic study of the phe- 
nomenon of IS in the spacetime of a single boson star, 
without a self-interaction potential for the non-minimally 
coupled scalar field. We choose to work directly in the 
physical Jordan frame where the physics is better un- 
derstood. In addition, we study some properties of the 
gravitational scalar waves such as their magnitude and 



frequency. All the numerical evolutions are performed 
in spherical symmetry using a 3+1 formalism of STT 
as presented in (24|. However, instead of evolving the 
geometry with the Arnowitt-Deser-Misner (ADM) equa- 
tions, we use a strongly hyperbolic version similar to the 
Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formula- 
tion [2a Hi| but adapted to the STT [27j]. Using this 3+1 
system the initial value problem of STT in the Jordan 
frame turns to be well-posed. Nevertheless, at the mo- 
ment we only have numerical evidence to show that the 
initial boundary value problem (IBVP) is also well-posed. 
The analysis of the continuum IBVP for this system will 
be presented elsewhere. 

This paper is organized as follows. Section |TT] intro- 
duces the STT and discusses briefly some properties as- 
sociated with the Jordan frame. The relevant 3+1 equa- 
tions of [24| are also presented. For completeness, and 
for the benefit of the reader, in Section llHl we discuss the 
heuristic analysis performed by Damour and Esposito- 
Farese in [l3| , which allows one to understand the scalar- 
ization phenomenon on simple grounds. Section [IV] con- 
tains our boson star model. In Section|V]we describe the 
scalar waves predicted in STT. Section I VII summarizes 
the setup used in the numerical code. We present the 
results of our numerical simulation in Section IVII1 and 
we conclude in Section IVIII1 

We also present in Appendix [A] the characteristic de- 
composition of the spherically symmetric equations used 
in this paper, which allow us to conclude that our system 
is strongly hyperbolic and, therefore, that the Cauchy 
problem is well-posed for the spherically symmetric case. 
In Appendix |B] we present some numerical evidence that 
indicates that corresponding IBVP is also well-posed. 

II. SCALAR TENSOR THEORIES OF GRAVITY 

A. Field equations 

The STT of gravitation are one of the simplest and 
most analyzed alternative theories of gravity. These al- 
ternative theories were introduced by Jordan during the 
decade of the fifties Q, and then reanalyzed by Brans 
and Dicke later Q. The general action for STT in the 
Jordan frame, where gravity is coupled non-minimally to 
a single scalar field <fi, is given by 

S[g ab ,<t>M = j {-^^ Q R -\^<t>? -V {<!>)} V^gd A x 

+ S ma , tt \g a b,ip] , (2.1) 

where ip represents all the matter fields, i.e. fields other 
than <j>, Go is the usual gravitational constant, F(<fi) is 
some NMC function to be specified later, and V{4>) repre- 
sents a potential for cj> (we use units such that c = 1). In 
fact, in all the numerical analysis considered here, we will 
not consider the potential V(<f>). However, for complete- 
ness it will be included in the field equations displayed 
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below. Notice that one can identify the "effective" grav- 
itational constant as the coefficient G c k(4>) = G /F(<p) 
that appears in the above action (63j |. 

From the above action, one finds the following field 
equations: 
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g a bR — 8TrGoT a b , 

f'R = V , 



(2.2) 
(2.3) 



where a prime indicates d$, □ := g ab V a^b is the stan- 
dard covariant d'Alambcrtian operator, and 



Tab 



Geff / / 

G 1 ab 



J"? 



Timatt 1 
ab J 



Tib ■= Va (/'V 6 0) - g ab V c (/'V c 0) , 



-(V0) 2 + y(0) 



(2.4) 
(2.5) 
(2.6) 



with T™ att the stress-energy tensor of all matter fields 
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Notice that equation (|2.2[) implies that the Ricci scalar 
can be expressed in terms of the trace of the energy- 
momentum tensor (|2.4p . Therefore, Eq. (|2 .3[) can be 
rewritten in the form 



□ c 
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with T mat t the trace of T^ att . On the other hand, the 
Bianchi identities directly imply 



V C T C 



0. 



(2-9) 



Nevertheless, the use of the field equations leads to the 
conservation of the energy-momentum tensor of the mat- 
ter alone 



v c-* matt — u : 



(2.10) 



which implies the fulfillment of the (weak) equivalence 
principle, i.e. test particles follow geodesies of the metric 

9ab- 



B. 3+1 decomposition 

In order to recast the previous field equations as a 
Cauchy initial value problem [2^, [29|, we first rewrite 
the four-dimensional metric in 3+1 form as 

ds 2 = -{a 2 - ft fa) dt 2 + 2 (3 t dx l dt + ~,,<l. <■■,!.<■■■ , (2.11) 

with a the lapse function, /3 Z the shift vector, and jij the 
3-metric induced on the spatial hypersurfaces. 



We perform a 3 + 1 decomposition of equations (|2.2[) 
and (|2.3p . using the normal timelike vector n a to the 
spacelike hypersurfaces £{, and the projection operator 
P a b = S a + n a rif,. In order to do so, we first define the 
first order variables: 



Qi := Di<f> = ff V fc < 



n := n a V a = 



dt ' 



(2.12) 
(2.13) 



where Di is the covariant derivative compatible with the 
3-metric 7^, and dj dt := dt — £p, with Cp the Lie deriva- 
tive along the shift vector. Notice that the relevant com- 
ponents of the quantities computed with P a b are the spa- 
tial ones. It is now straightforward to show that Qi and 
IT evolve according to |24( , 
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We define the energy density p := n a n b T a b, momen- 
tum density J a := —P b a n c T 0C and a stress tensor 
Safe := P c a P d bT cd . From (El we find 
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Using now Eqs. (|2.7|) and 
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where we have defined Q 2 := Q l Qi- 

Notice that the 3 + 1 decomposition of (|2.2j) are just 
the usual ADM equations given by 



dt 
dKjj 
dt 



-2aKi, 



(2.22) 



-ViV 3 -Q + a [Rij + KKij - 2K ll K l ] ] 
AnGoa (S - p) - 2S t] ] , (2.23) 



where Rij is the 3-dimensional Ricci tensor associated 
with the spatial metric 7y , and the effective matter terms 
are given by Eqs. (|2.19|1 - (|2.21|1 . The Hamiltonian and 
momentum constraints take the form 

H := i (R+ K 2 - K i:j K ij ) - 8nG p = ,(2.24) 



M % := Di (K a - Y K) ~ 8ttG J 1 = 



(2.25) 



Formally, one should also consider the constraint 
D[iQj] — 0, which corresponds to the integr ability con- 
dition dfj(f> = d 2 ^. The above system of evolution equa- 
tions has to be completed with appropriate evolution 
equations for the gauge variables. This issue is considered 
below. 



C. Gauge choice 

To obtain a closed evolution system, one has to im- 
pose gauge conditions for the time variable t and for 
the spatial coordinates x 1 . Following [27]], we will con- 
sider a modified Bona-Masso (MBM) time slicing con- 
dition (30|. In local coordinates adapted to the 3 + 1 
foliation x a = (t^x 1 ), this slicing condition is given by 



da 
~dt 



K- 



e /' 



fBAt(a) I 



n 



(2.26) 



with /bm(o) > the usual Bona-Masso gauge function 
and 8 a free parameter. 

The specific choices G = /bm = 1 correspond to a 
modified harmonic slicing condition (termed "pseudo- 
harmonic" in [3, HtJ ) , which was specially useful for the 
second order hyperbolicity analysis performed in [24| . On 
the other hand, with = one recovers the usual Bona- 
Masso (BM) slicing condition. However, it has been 
shown that taking = does not result in a stro ngly 
hyperbolic formulation of STT in the Jordan frame 27 1 . 
For this reason, in all the simulations presented here we 
have used the pseudo-harmonic slice with = 1. 

Concerning the propagation of the spatial coordinates, 
we will consider the shift vector as an a priori known 
function of the coordinates. In particular, in all our evo- 
lutions it is set to zero. However, in the future, it would 
be interesting to investigate some "live" shift conditions 
and their effects in phenomena involving STT. 



III. SCALARIZATION 

The STT of gravity induce strong field effects which, 
for instance, produce important deviations from GR in 
stellar models. As mentioned before, one such effect is 
the SS phenomenon which is similar to the spontaneous 
magnetization in ferromagnetic materials at low temper- 
atures. In the following, we will use this analogy in order 
to understand the SS phenomenon. 

When a ferromagnet is exposed to an external mag- 
netic field, the individual spins of its constituents align 
with the field, giving rise to a permanent magnetization 
which remains even after the external field is switched 
off. Moreover, the ferromagnets have the property that, 
below the Curie temperature, a magnetization appears 
"spontaneously" even in the absence of an external mag- 
netic field. In the STT, on the other hand, a nontrivial 
configuration of a scalar field may spontaneously appear 
during the evolution of a compact object in the absence 
of external sources, i.e. without a potential V{<p). One 
can then identify the external magnetic field with a back- 
ground (cosmological) scalar field and the temperature 
with the inverse of the central energy-density /?™ att of the 
matter content or, equivalently, with the inverse of total 
baryon mass (in the case of neutron stars). The roll of the 
magnetization is played by a new global quantity called 
the scalar charge Q sca i, which will be defined below and 
which corresponds to the coefficient of <p(r) ~ Q S cal/f i n 
the asymptotic region. This means that beyond a cer- 
tain critical density or critical baryon mass, the transi- 
tion to the spontaneous scalarization ensues. In this case 
dQscai/dp™^ is infinite at the critical energy-density. 
This transition can be smoothed by the presence of a 
non-zero background scalar field <f>Q. 

In practice, when the phenomenon is analyzed in static 
configurations, the value 4>o is usually fixed by a shooting 
method [H, [3l[ . It is important to emphasize the rele- 
vance of not adding sources to the scalar field equation. 
Indeed, a non-zero potential can make the scalar field de- 
crease so fast that the scalar charge might in fact vanish 
at infinity. 

An interesting analytical toy model to understand this 
phenomenon is the following Q. Consider a static and 
spherically symmetric compact object represented by 
an incompressible fluid (constant energy-density), whose 
profile density is given by a step function. Moreover, 
assume that the function F((j>) in the action for STT 
is just a quadratic function of, F(</>) = 1 + 8irGo(;<f) 2 , 
and the potential V((j>) vanishes. Finally, assume R « 
-8 7r G Tmatt ~ 8 7r G p matt (as in GR). This implies, 
for instance, that the generalized Klein-Gordon Eq. (|2.3|) 
is linear in </>. Furthermore, by neglecting all the grav- 
itational effects in this equation, one ends-up with a 
Hclmholtz-like equation of the form A<f> + k 2 (j> = 0, where 
the mass-like term k 2 = 8 ir Go £ p ma tt depends on both 
the constant energy-density of the matter and the NMC 
constant £. Notice that in this simplified model, the mass 
term vanishes beyond the surface of the compact object. 
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Now, for £ > the interior regular spherically symmetric 
solution of the above Helmholtz equation is given by 



<knt(r) = 4>, 



sin(fc r) 
kr 



(3-1) 



where <f> c is the scalar field at the origin r = 0. The 
exterior solution is 



C 

4>ext{r) = h 0o 

r 



(3.2) 



where C is an integration constant and O is the asymp- 
totic value of the scalar field. 

When both solutions are matched continuously at the 
surface of the object r = 1Z [p ma tt(f > Tty = 0], it 
turns out that <f> c = <f>o/ cos(klZ). The explicit form for 
the constant C is not relevant for the analysis (but one 
finds C oc 4> C 1Z). Note that if <f>o = 0, automatically 
<j) c = = C and, therefore, <j>(r) = 0. In this case, 
there is no scalarization. On the other hand, a different 
situation can happen if cos(fc7?.) vanishes as well when 
4>o —> 0. This can occur when k = n/(2TZ). In this case, 
the solution is given by 



4>mt(r) = 

0ext(»") = 



e sin(r) 



(3.3) 
(3.4) 



where e is a constant related to the scalar charge Q SC ai 
whose numerical value depends on the details of the 
model, and f = rn/(21Z). The above simplified anal- 
ysis agrees qualitatively with the full numerical study 
H, [HI, [UJ]. Note that if £ < 0, the interior solution is 
<knt(r) = c sinh(|fc| r)/(\k\ r) and 4> c = <f> / cosh(|fc| 11). 
In this case, the scenario is completely different since the 
function cosh(|fc| TV) never vanishes, and so when 4>q — > 
then automatically <p c — ¥ and then the scalar field van- 
ishes everywhere (no scalarization ensues). It is some- 
how remarkable that the spontaneous scalarization phe- 
nomenon is associated with a decreasing effective gravi- 
tational constant {i.e. G e s < Go) [l8j ]. 

Another way to understand the existence of these kind 
of scalarized configurations is to notice that the presence 
of a non-trivial scalar field <f>(r) (within a class of STT) 
causes the total energy of the stationary configuration 
to decrease relative to the case where <fi(r) = for a 
fixed baryon mass [HI, [l8[ . This can also be understood 
on Newtonian grounds by a suitably redefinition of the 
kind of energy that has to be minimized when dealing 
with a theory where an effective gravitational "constant" 
may vary [l8|. The energetic analysis shows that for 
large compactness, the energetically preferred stationary 
configurations are those with a non-trivial scalar field. 
Again, in the ferromagnetic analogy, one appreciates that 
in the Landau ansatz, the free energy of the ferromagnet 
becomes lower in the presence of magnetization than the 
energy in the absence of it when the temperature is below 
the Curie point. This occurs since below that tempera- 
ture the free energy develops a global minimum and a 



local maximum (like a Mexican hat potential). The local 
maximum of the free energy is located at zero magnetiza- 
tion while the minimum corresponds to a non-zero mag- 
netization. Recently, it was also been argued that the 
SS phenomenon can be traced back to the quantum fluc- 
tuations of the vacuum state associated with the scalar 
field [Ull. 

The scalar charge which characterizes the scalarized 
configuration is defined as 



Q 



seal 



:= — lim 
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s a S7 a (j)ds 



(3.5) 



where s a is the unit outward normal to a topological 2- 
sphere S embedded in S t , and r is a radial coordinate 
that provides the area of S asymptotically. As it was re- 
marked in the introduction, when the asymptotic value of 
the scalar field, 4>q, is not demanded to vanish but is only 
accommodated to satisfy the Solar System bounds, then 
the scalarization process is induced by such background 
(cosmological) field. In such situation the transition from 
a small scalar charge to a large one (which depends on 
the compactness of the object) is smoothed out and the 
derivative dQ sca i / d p natt is always finite. In this paper, 
we will be concerned with this latter situation only, but 
as long as (f>o is small, the difference between the two type 
of scalarizations is just a matter of principle. Neverthe- 
less, the important point for making such a distinction is 
that while in the SS case there are no emission of scalar 
gravitational waves (when i 7 "((/))0 o= o = 0, which is the 
case for the quadratic function F{_4>) considered above), 
for the induced case where 4>o 7^ 0, one can have a small 
but non-zero amplitude for the scalar waves (c.f. Sec. lV)) . 

Another important qualitative mathematical aspect 
that distinguish both types of scalarizations for a NMC 
like the quadratic one is the following. If one consid- 
ers Eq. (|2.3p in absence of a potential, it turns out that 
<f> = is always a stationary solution of the equation. 
This implies then that 0o = 0. Therefore, in order to 
trigger the transition to a SS case an explicit scalar-field 
perturbation is required. An analysis of this sort was 
performed by us in [2p| . However, if one considers ini- 
tially a trivial (but non-zero) scalar-field configuration 
<p — 4>q = const., then this is not a stationary solution 
of Eq. (|2.3[) . In such a case a fortiori the scalar field 
will evolve in time without the need of any perturba- 
tion. How much of the initial energy of the star will then 
be transformed into scalar energy and scalar radiation 
leaving behind a highly non trivial stationary scalar field 
configuration will depend precisely on the compactness 
of the object. Higher compact objects will radiate more 
energy in the form of scalar radiation than lower compact 
ones. Therefore, higher compact objects are expected to 
end up in a stationary scalarized state with energy lower 
than the initial one, the difference being radiated away 
in scalar-field form. Clearly, in order to analyze in de- 
tail the transition towards a scalarized state from a state 
with a trivial (non-zero) scalar field and its correspond- 
ing emission of scalar radiation, a dynamical evolution is 
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required. This is the aim of this paper. 



IV. BOSON STARS 

Boson stars are equilibrium configurations of a self- 
gravitating (condensate) complex scalar field. Their "hy- 
drostatic" equilibrium is maintained by the intrinsic ef- 
fective pressure of the boson field due to the uncertainty 
principle (for a review see Ref. [HI, [Hf| ) , rather than the 
Pauli exclusion principle that acts, for instance, in neu- 
tron stars. Classically, one can interpret the equilibrium 
as a consequence of an effective pressure associated with 
the boson field which depends on its gradients and po- 
tential. Since the energy-density and pressure are pa- 
rameterized in a certain way by the boson field itself, the 
relation between them provides a non-trivial EOS for this 
kind of matter. 

Boson stars are also interpreted as macroscopic boson 
quantum states whose associated physical particles are 
formed by the excitations around the vacuum expecta- 
tion value of the scalar field. The theoretical existence 
of such objects were proven first by Kaup [HI, and lat- 
ter by Ruffini and Bonazzola [37j . for the ground state 
solutions of a free boson field. Using the uncertainty 
principle and the definition of the Schwarzschild radius, 
it can be shown that the boson stars considered there 
have an effective radius R e ff ~ h/rribC and a maximum 
mass of 
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2 G m b 



0.5 MlJm b , 



(4.1) 



where Mpi is the Planck mass and m b is the mass associ- 
ated with bosons (for clarity we have restored the speed 
of light c). Numerical results show that, in fact, this limit 
is M max » (2/V) Mpjmb- Therefore, the resulting sizes 
and masses of boson stars would be so small as to be 
astrophysically inconsequential. In a more recent paper 
[HI , it was shown that a self- interacting boson field (with 
an interaction of the form ~ A</> 4 ) can give rise to stable 
boson stars with 



M max ~ X 1 ' 2 M 2 Pl /m b ~ GcV 2 X 1/2 M Q /m 



(4.2) 



which is comparable with the Chandrasckhar mass for 
fermion stars [3|| . This result was then extended to the 
so-called soliton stars [4(| Hl| . 

Boson star models have been constructed in the pas t 
within the framework of STT (see e.g. 0, H3, 1H- 
and references therein) , although only Whinnett had 
shown that the phenomenon of spontaneous scalarization 
occurs in these objects with the inclusion of a quartic 
self- interaction potential [l9j|. Recently, we have found 
that self-interactions are not in fact necessary in order 
to produce scalarization [2(|. On the other hand, Tor- 
res has shown in [42j ] that, for parameters and boundary 
conditions respecting the weak-field limits and the nu- 
cleosynthesis bounds, the masses of boson stars in the 



STT framework are comparable with the ones in GR 
(for stars in the ground state). Comer (43[ confirmed 
the same trend for the case of boson stars in "excited" 
states. Equilibrium and stability properties for these 
stars in STT for different cosmic ages have been analyzed 
in [MSI]. 

An important aspect of analyzing boson stars in the 
framework of STT is that the transition to a scalarized 
state might be accompanied by the emission of (spin- 
0) scalar gravitational waves (like in neutron stars) . It is 
possible that such kind of waves might be detected in the 
future if fundamental scalar fields do exist in nature 12 ill . 



A. The Model 

Boson stars are described by the Lagrangian density of 
a complex scalar field 



i 



<7 Qb V a ^V & ^-^(|V| 2 ): 



(4.3) 



where tp is the scalar field, ip* its complex conjugate, 
\ip | 2 = ip tjj*, and ^(l^l 2 ) is a potential depending just 
on the norm. Variation of the above Lagrangian with 
respect to ip leads to the KG equation 



Dip = 2tp 



(4.4) 



It is straightforward to show that if the scalar field 
ip is real, the KG equation (|4.4p takes the usual form 
Uip = d^pVrt,. On the other hand, variation of (|4.3j) with 
respect to the metric g ab leads to the energy-momentum 
tensor 



rTimatt 
1 ab 



1 r 
2 



(V a ^*)(V h ^) + (v 6 v*)(VaV0 



~9ab 



\m\ 2 



(4.5) 



We will consider only the free field case where the poten- 
tial is given by 



1 



(4.6) 



with m a parameter that can be consider as the "bare 
mass" of the theory (although it has units of inverse 
length). It is possible to include more general terms 
in (|4.6p . such as X\ip\ A which corresponds to a self- 
interaction term [38j . 

The Lagrangian (|4.3[) is invariant with respect to a 
global phase transformation ip — > e lq ip (with q a real 
constant). The Noether theorem then implies the local 
conservation of the boson number V ' a J a = 0, where the 
number-density current is given by 



(4.7) 
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This means that the total boson number, 
Af = - 



(4.8) 



is conserved (here 7 is the determinant of the 3-metric 
7ij). The total boson mass can be defined by 



rribAf , 



(4.9) 



where m& := 2 ir hm/c is the mass of single bosons (again 
we have restored the speed of light c) . 

A 3 + 1 decomposition of the above energy-momentum 
tensor gives rise to the following matter variables 



Pmatt^i (|rM 2 + |Q V ,| 2 )+^(|Vf), 



°matt 



J? 



-(|Qvl 2 -|n^| 2 ) + Wl 2 ) 



(4.10) 

(4.11) 
(4.12) 



where we have defined, in analogy with (|2.12[) and (|2.13l) , 
the variables 



Qt 

U 4- 



1 dip 
a dt 



a v 



(4.13) 
(4.14) 



According to this, one can rewrite the KG equation (|4.4[) 
as a first order PDE system like (|2. 14112. 151) . 

Before finishing this Section we must emphasize the 
fact that boson stars in STT involve two distinct scalar- 
fields, the real- valued non-minimally coupled field </>, and 
the a complex- valued boson field tp. 



V. SCALAR GRAVITATIONAL WAVES 

In this section we will show that the presence of 
the scalar field <fi can induce the propagation of scalar 
(monopolar) gravitational waves. Similar analysis have 
been presented before in [U |48| . Using the weak- field 
approximation, Wagoner has analyzed the properties of 
the source and its radiation in the Einstein frame j47j . 
On the other hand, Harada et al. have analyzed, in the 
Brans-Dicke theory, the emission of the scalar gravita- 
tional radiation in spherical dust collapse [48[ . Moreover 
from the detection point of view, a detailed study has 
been performed in |2l| . 

Let us start by considering a linear perturbation of a 
fiat background metric rj a b and a background scalar field 
4>o such that 



9ab ~ Vab +ejab, 



(5.1) 
(5.2) 



where e <C 1 (do not confuse this e with the one in 
Eqs. (|3.3p and (|3.4p ). Thus, according to the above ap- 
proximation, we have 



Tab ~ e T a b , 



(5.3) 
(5.4) 



where the subindex indicates quantities at zero order. 
Notice that one must in fact have Ttf b = in order to be 
consistent with first order perturbations of a flat back- 
ground (49l | . 

In the generalized Lorentz gauge d a i a b = 0, where i a b 
is defined as 

lab ■= lab -^f] a b(l + 2 4> j^j , (5.5) 



the field equations flUD and flU} ( for V = 0) become [49| 

(5.6) 
(5.7) 



□ f! 7o6 = -167r^2r* t 1 



□ , 



F' ~ 
4-?rC — T mat t , 
^0 



where D v is the D Alambertian operator in the flat back- 
ground metric, and the constant £ is defined as 



3(F, 



(5.8) 



Notice that the flat metric is used to raise and lower 
indices of first order tensorial quantities. 

We will now perform the analysis of scalar gravita- 
tional waves in a spherical and vacuum spacetime, which 
is enough for the purpose of the following discussion. 
In that case one can neglect the tensor modes and as- 
sume that 7 a b = (note that in spherical symmetry, 
quadrupole radiation is absent). Hence, according to 
equation (|5.5[) we obtain 



(5.9) 



lab = -~4>f]ab "pr ■ 
^0 

Therefore, the whole metric at linear order reads 
9ab ~ (1 + $ ) Vab , 



(5.10) 



where $ := —4>Fq/Fq and the factor e has been re- 
absorbed in (j> (note that <f> *C <j>o (*- e - ^ "C 1))- The 
physical metric then turns out to be conformally flat, 
g a b ~ f2 2 T] a b, with the conformal factor O 2 := 1 + $. 

Gravitational waves are directly related with the Rie- 
mann tensor. Here we compute that tensor using the 
well-known relationship between two Riemann tensors as- 
sociated with two conformal metrics (see e.g. Eq. (D.7) of 
Wald's [5CiJ]). In this case, the Riemann tensor associated 
with rjab vanishes and, at linear order, we have 

RLab/ = S d [a V b] V c $ - lf e T} c{a V 6] V e $ , (5.11) 
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where the subindcx L indicates that this is valid only at 
the linear approximation. According to (|5.1ip , we obtain 
the following components, which are directly related with 
the relative (tidal) acceleration between two particles in 
geodesic motion, 



Rl 



titj 



(5.12) 



VI. NUMERICAL SETUP 

In this section we present the numerical ingredients 
that have been used in order to study the scalarization 
transition and the emission of scalar radiation in the bo- 
son star context. 



Assuming now that <f> = <j)(t, r) is a spherically sym- 
metric perturbation for the scalar field, we obtain 



Rl 



-6ijd? t <$> + SiSjd 2 r <f> + - (Sij - SiSj) d r <5> 



. (5J3) 

where s, = SijX 3 /r is a unit vector in the radial direction 
of propagation. 

On the other hand, in vacuum and in spherical sym- 
metry, the wave equation (|5.7[) implies 



d 2 J = d 2 J + - d 7 . 



(5.14) 



Plugging (pTT4l into (f5TT5|) yields 



titj 



- (Sij - 3siSj) 



(5.15) 



For outgoing radiation r) = ty(t — r)/r, thus the sec- 
ond term in (|5.15|) . which involves the spatial derivative, 
will be very small with respect to the first one in the 
"wave zone" and so we can neglect it. Finally, we obtain 
the following expression 



Rl 



titj 



IJ u tt 



(5.16) 



n 



where we have introduced the transverse projector 
J-ij = Sij — SiSj in the orthogonal directions to Sj and 
neglected terms 0(l/r 2 ). In the above expression, we 
have considered only the massless case. When the mass 
term is included, a longitudinal contribution appears in 
Eq. (I5U61) (c.f. Ref. [III). 

According to (|5.9p . in the minimal coupling case, where 
0, or when one takes the value </>o = asymptoti- 
cally (like in the quadratic coupling F = 1 + 87r Go£</> 2 ), 
the scalar-gravitational waves are absent. However, tak- 
ing for 4>q the maximal value allowed by the Solar System 
experiments, which constrain wbd > 4 x 10 4 , scalar grav- 
itational waves are expected to develop in dynamical sit- 
uations where the NMC scalar-field is non-zero initially, 
like in the IS phenomenon. We expect the amplitude of 
the 1 /r-contribution of the gravitational radiation to be 
of the following order 



\Rl 



utj | 



Qscal W 2 

2D 



I -L; 



\K/Fo\ 



(5.17) 



where uj is the frequency of the scalar wave, Q sca i is the 
scalar charge (which has units of mass) and D the dis- 
tance to the source. 



A. Formulation 

We perform a numerical evolution of the equations of 
Sec. Ill Bl However, since the ADM equations in GR are 
only weakly hyperbolic fo"H , we will use a formulation of 
the BSSN type, which has been particularly robust in the 
numerical evolution of both vacuum and matter space- 
times in GR [25|, [26| . We adopt the particular version of 
the BSSN formulation presented in 52| that is specially 
adapted to spherical symmetry. In Appendix [A] we dis- 
cuss the characteristic decomposition of this formulation 
for the case of an STT. 



B. Initial Data 

We consider a single boson star in stationary equilib- 
rium. In such a configuration the spacetime metric is 
time independent, and the scalar field ip(t,r) oscillates 
in time with a fixed frequency u: 



ip(t, r) = *(r) e l 



(6.1) 



In order to find the initial data one must substitute (|6.1[) 
into the KG Eq. (|4.4p . We now need to solve Eqs. (|2.2p 
and (|4.4p in order to obtain the frequency ui, the func- 
tion ^(r), and the metric coefficients such that, for a 
given amplitude of the scalar field at the origin, ^(0), 
the resulting spacetime is static. Following |36l. [53l [54|. 
we solve this problem in the polar-areal gauge, where the 
line element takes the form 



ds 2 = 



\r)dt 2 + A(r)dr 2 + r 2 dtf 



(6.2) 



where dVL 2 = d6 2 + sin(#) 2 dej) 2 is the usual solid angle 
element. The field equations then reduce to: 



d r A = A 



I- A 



Auj 2 ^ 2 
Airr { 2AV 4 , + ^ 



d r a 



A- 1 d r A 
r + ~2A 



%-ktAV, 



-Qv 
A^ 



2 d r a d r A 
r a A 



dibVib T 

or 



(6.3) 

(6.4) 
(6.5) 



(6.6) 
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where the potential is given by the Eq. (|4.6[) . There 
are several solutions of the above system, depending on 
the value of the different variables at the origin and their 
asymptotic behavior. In order guarantee regularity at 
the origin we impose the following boundary conditions 

A(0) = 1, 
d r a\ r= o = , 
*(0) = * c , 
Qv(0) = 0, 

Also, for the spacetime to be asymptotically flat we must 
ask for 



A\ XimT _ 

Ot\\ira r _ 



= I- 
= I- 
= 0. 



(6.11) 
(6.12) 
(6.13) 



One can consider this problem as an eigenvalue prob- 
lem for \P(r) in the sense that, for a given value of the 
scalar field at the origin the above system only ad- 
mits solutions for a discrete set of frequencies u>. We are 
interested in the ground state of the boson stars, which 
corresponds to the configuration with no nodes in *5?(r). 

Numerical solutions of Eqs. (j6.3H6.6p satisfying the 
boundary conditions (|6 .7116.1 31 . are obtained for a given 
value of V&c by integrating the equations from the origin 
outwards using a shooting method |3lj . This is similar 
to the case of neutron stars, where the integration is also 
performed outwards from the origin by giving a value of 
the energy-density at the center of the star. Fixing the 
value of the energy density at the center of the star is 
equivalent to choosing a particular stationary configura- 
tion, which in turn determines the ADM mass as well as 
the total number of particles (total baryon number in the 
case of neutron stars and total boson number in the case 
of boson stars). 

Moreover, just like in neutron stars, one can construct 
a whole family of stationary configurations for different 
energy-densities at the origin. It can be shown that there 
is a specific value of the central density which maximizes 
the ADM mass. That point indicates the threshold to 
configurations that are unstable to gravitational collapse 
under small perturbations. Figure [1] shows the mass pro- 
file for a single boson star without self-interaction. The 
maximum mass configuration, which corresponds to a 
central value of the boson field given by er(0) ~ 0.272 
[where a(r) := y/4 n G *&(r) is a dimensionlcss boson 
field; c.f. Eq. (|7.2[) below], separates the space of config- 
urations into two branches: the stable "S-branch" , and 
the unstable "U-branch)". When perturbed, configura- 
tions on the U-branch can either collapse to form black 
holes or disperse away depending on value of the binding 
energy (55j . 

In the case of STT, there is another point of instabil- 
ity towards spontaneous scalarization for a mass which is 
not the maximum one. In neutron stars there is a critical 
baryon mass below the maximum that marks the onset 



/ 

i 



(6.7) 


E 


(6.8) 


<0.4 ; 


(6.9) 


1 

:i 


6.10) 


j 

0.3 +■ 



Maximum mass 



U-branch 



S-branch 



0.25 0.5 0.75 



1.25 
o(0) 



FIG. 1: ADM mass profile for a single boson star without self- 
interaction in GR, as a function of the central value of the 
boson field <r(0). The maximum mass configuration, which 
corresponds to <r(0) ~ 0.272, separates the stable configura- 
tions (S-branch) from the unstable ones (U-branch). 



of this instability. We expect that it is the total boson 
number that indicates the point of instability towards 
spontaneous scalarization. However, a priori this value 
is difficult to know unless one constructs a whole fam- 
ily of stationary configurations and looks for solutions 
with a non-trivial value of the non-minimal scalar field 
cj). At this point is perhaps important to remark that a 
maximum mass model within STT usually corresponds 
to a model with a non-trivial scalar field. That is, it 
corresponds to a star that has already undergone a spon- 
taneous scalarization process in which some of the scalar 
field has been radiated away and the star has reaches 
a stationary configuration. This maximum mass star is 
then unstable and collapses to a black hole under a small 
perturbation (see Section |VII[) . 



VII. NUMERICAL RESULTS 

We performed numerical simulations of several single 
boson star configurations. We used the code described 
in [|52[, which consists of a spherical reduction of the 
BSSN formulation, coupled to the STT given by the La- 
grangian (|2.1[) . We have taken the non-minimal coupling 
function (12.71) to be of the form 



F{(j>) = 8tt/(0) = 1 + SttGo^ 2 



(7.1) 



with £ a positive constant. In order to have a notation 
consistent with previous studies about boson stars (42l . 
l5a. l56fl , we define the dimensionlcss boson field 



(7.2) 



All runs have been performed using a method of lines 
with a fourth order Runge-Kutta integration in time, 
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and fourth order centered differences in space. Con- 
straint preserving boundary have been implemented us- 
ing the algorithm described in Appendix El Our typi- 
cal simulations use a grid spacing of Ar = 0.048, and 
we take At = Ar/2 in order to be sure that we satisfy 
the Courant-Friedrichs-Levy stability condition [HT| ■ We 
have also performed some simulations with grid spacings 
of Ar = 0.024 and Ar = 0.012, in which case the vio- 
lations on the Hamiltonian constraint are dominated by 
spurious reflections from the boundary which are, in the 
worst case, of order 10~ 9 . In all the runs presented here 
the outer boundaries are located at r out = 240. 



A. Spontaneous Scalarization 

The phenomenon of scalarization in compact objects 
and the emission of the monopolar gravitational waves 
depends strongly on the asymptotic value of the non- 
minimally coupled scalar field cf>. If the asymptotic value 
of <j> is zero, which corresponds to spontaneous scalariza- 
tion, then F' tends to zero asymptotically which means 
that, according to (|5.17[) . there are no monopolar grav- 
itational waves. Nevertheless, it has been shown that 
in this case the evolution of stable boson stars on STT 
reaches a final state where the scalarization ensues poj . 
Boson star configurations below a critical central value of 
the boson field ^ !lt are stable with respect to Gaussian 
perturbations on the scalar field and do not lead to a SS 
transition: The scalar field cf> is just radiated away during 
the evolution. On the other hand, configurations above 
that critical value are unstable with respect to perturba- 
tions and undergo a transition to a scalarized state with 
a non-trivial scalar field <fi and a non-zero scalar charge. 



parameter cjbd can be written as 

_ 1 + 87tG £$ 



J_ 

f,2 



32ttGoC 



2^ 



(7.3) 



where £ is the NMC constant associated with the class 
of STT considered here, it turns out that 



1 



v/8< (4w BD £- i; 



(7.4) 



For a given value of £, one can then enforce the constraint 
wjjd ^ 4.3 x 10 4 imposed by the Cassini probe [13, and 
use the above expression to fix the value for (f> that sat- 
isfies the observational bounds. In all the simulations 
discussed below we have chosen </> in such a was as to 
saturate the Cassine bounds, that is, the one correspond- 
ing to ojbd = 4.3 x 10 4 (notice that </>o will then depend 
on the parameter £). 
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B. Induced scalarization 

In the IS phenomenon, an initially non-zero and uni- 
form NMC scalar field evolves naturally without the need 
of any perturbation. The reason for this is that for a 
non-zero background scalar field <fi(r) = (f>o one finds 
F' {4>)\<l>=<i>o 7^ 0. The presence of the term f'R in the 
Klein-Gordon equation (|2.3[) then implies that this ini- 
tial data will not be a stationary solution, and the field 
will therefore evolve away from its initial configuration 
without the need for an external perturbation. This is 
in contrast with the SS case for which we have initially 
4>{r) — 4>q — 0, which is indeed a solution of Eq. (|2.3I) . 
so an explicit perturbation is required to trigger the SS 
phenomenon. 

Since in principle any background scalar field </>o in 
STT might disturb the constraints imposed by the Solar 
System experiments, the value of <f>o must be chosen such 
that the corresponding Brans-Dickc parameter lobb sat- 
isfies the observational bounds. Using the fact that the 



FIG. 2: Snapshots of the evolution of the NMC scalar field 
for the case when <r(0) ~ 0.266 initially. For this simulation 
we have taken £ = 1, and the NMC scalar field is initially set 
to 4>{r) — 4>o = —4.8 x 1CP 4 , which corresponds to the back- 
ground value that saturates the lower bound for lobd- The 
system evolves until it reaches a quasi-stationary configura- 
tion with a non trivial scalar field tj>. This final configuration 
is what we refer to as "induced scalarization". Notice that 
asymptotically <j> preserves its initial value 0o- 

Figure [2] shows snapshots of the evolution of the NMC 
field <p for a boson star with an initial central density 
cr(0) ~ 0.266, which is at the threshold of the unstable 
configurations. For this simulation the NMC parameter 
was taken to be £ = 1, and the NMC scalar field is ini- 
tially set to 4>{r) = <fio = —4.8 x 10~ 4 , which corresponds 
to the background value that saturates the lower bound 
for ojbd- In this case a quasi-stationary configuration 
with a non-trivial scalar field 4> is reached at the end of 
the simulation. 

Figure [3] displays the initial and final states of the ra- 
dial metric component A, the lapse function a, the NMC 
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FIG. 3: Initial and final states of the radial metric component 
A (top left), lapse function ot (top right), the NMC field 
(bottom left), and the norm of the complex bosonic field ^ 
(bottom right), for an initial central density er(0) ~ 0.106 and 
£ = 500. Just as in Fig. [2j the system undergoes scalarization 
and reaches a quasi-stationary configuration with a non-trivial 
NMC scalar field at the end of the evolution. 



scalar field 0, and the norm of the boson field |^|, for a 
different simulation with £ = 500 and <r(0) = 0.106 ini- 
tially and 0o = —9.6 x 10~ 7 . Notice that, although the 
parameter £ is much larger than in the previous case, the 
IS only seems to affect slightly the geometry of the space- 
time (confront, for example, the initial and final values 
of the lapse ct(0) in Fig. [3]). In order to make this state- 
ment more precise, we first define an approximate size of 
our boson star i? s t ar as the radius of the sphere which 
contains the 95% of the integrated mass. The compact- 
ness of the star is then defined as MMs/-Rstar, with Mms 
the Misner-Sharp mass of the system (see below). For 
the simulation of Figure [3] we find that the compactness 
of the star is reduced by only ~ 3% with respect to its 
initial value. 

The Misner-Sharp mass Mms used to measure the mass 
of the system is defined as follows: Let r a be the areal 
radial coordinate, we first define the "Misner-Sharp mass 
function" toms (fa) if* terms of the radial metric compo- 
nent gr a r a ( r a) aS [58[ 



m M s(r a ) ■= y fl 



9r a r a (Tcl) 



(7.5) 



with r a the areal radius. One can show that, 
for asymptotically-flat spherically-symmetric spacetimes, 
the Misner-Sharp mass function coincides with the ADM 
mass as long as we are in vacuum. That is, once we 
are in a region outside of all the sources we should find 
that mMs(r ) reaches a constant value Mms such that 
M M s = M A dm- 

Several considerations arc now in order. First, even 
though our initial data is in the areal gauge, during the 



evolution this is no longer the case, so that one must 
transform back from the radial coordinate r used in the 
simulation to the areal radius r a in order to calculate the 
mass function. This is simple to do and we will not go 
into the details here. More important, however, is the 
fact that for boson stars we are never actually in vacuum 
since the bosonic field 'F extends all the way to infin- 
ity. However, one finds that <F decays exponentially, so 
that in practice for stars that are not scalarizcd we very 
rapidly reach a region that for all practical purposes is 
indeed vacuum. For scalarized stars, on the other hand, 
we have to be more careful since the NMC field de- 
cays more slowly (typically as ~ and in the IS case 
it in fact reaches a small non-zero asymptotic value 0o- 
This means that the Misner-Sharp mass function never 
quite reaches the constant "vacuum" value. Because of 
this, the values of Mms wc report actually correspond to 
the value of the mass function toms(?") evaluated at the 
boundary of the numerical grid. Strictly speaking this 
is not the actual mass of the system, but it is a good 
enough approximation for our purposes. 

As a final comment one should also mention the fact 
that the Misner-Sharp mass Mms only coincides with the 
ADM mass for static spacetimes. The spacetimes con- 
sidered here are dynamic, and radiate energy during the 
scalarization process. When we talk about the initial and 
final mass we mean in fact the mass of the star, and not 
the "true" ADM mass which takes into account contribu- 
tions all the way to infinity (and is therefore constant). 

We will now try to understand the emergence of the 
scalarization phenomenon on energetic grounds, where 
stationary scalarized configurations turn out to be ener- 
getically preferred over unscalarized ones. For instance, 
it has been shown in the context of neutron stars [H, [H| 
that, beyond some critical baryon-mass, a stationary con- 
figuration with $1^0 which maximizes the fractional 
binding-energy of the system, is energetically more favor- 
able than the corresponding configuration at the same 
baryon mass with 0=0 (the GR case). This critical 
point depends on the details of the model, such as the 
value of £, the equation of state, etc. 

Figure [4] shows the results of this energetic analysis 
for our scalarized boson stars. In the Figure we plot the 
fractional binding-energy (Mdos/Mms) — 1 as a function 
of the bosonic mass Mb os for several values of £. Here 
Mms is calculated as described above, while Mb os is the 
mass contribution from the bosonic field alone obtained 
from Eq. (|4.9[) . From the Figure one can see that, for a 
given value of M^ os and £, the scalarizcd configuration is 
energetically preferred when compared with an ordinary 
boson star with 0=0 (the GR case). It is important to 
mention here that, strictly speaking, this energetic anal- 
ysis would only be valid for the SS case, where both the 
scalarized and unscalarized configurations have a vanish- 
ing asymptotic value of the NMC field 0. Nevertheless, 
for small asymptotic values 0o one finds that bulk quan- 
tities are essentially the same in the IS and SS cases, so 
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one can still use the energetic argument. 
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FIG. 4: Fractional binding energy (Mh os /MMs) — 1 as a func- 
tion of the total boson mass M^ OB , for both stationary scalar- 
ized and unscalarized (GR case) configurations. Configura- 
tions with larger binding energy are energetically preferred. 
The points correspond to the actual quasi-stationary final 
configurations resulting from our simulations, while the lines 
are a simple fit obtained with polynomial regression. 

To corroborate the above analysis, in Figure [5] we plot 
the final value of the integrated scalar charge Q sca i, and 
the central absolute value of the NMC scalar field <fi 
(which always corresponds to its maximum), for several 
sets of boson star configurations with different values of 
the parameter £ as a function of the initial central density 
of the star <r(0). We have chosen initial configurations 
with central density in the interval [0.028,0.265], which 
correspond to boson stars which in GR are stable against 
gravitational collapse (i.e. the S-branch configurations 
of Fig. [1]). For these configurations we have always as- 
sumed cjbd = 4.0 x 10 4 , which fixes the asymptotic value 
of the NMC field <f>Q. As expected, these stars evolve to a 
quasi-stationary final state with a non-trivial scalar field 
4>. We can conclude that the scalarization phenomenon 
depends strongly on the value of the parameter £. For all 
configurations that we have analyzed we find that the fi- 
nal value of the scalar charge, as well as the central value 
of (j>, reach a maximum for £ ~ 50. Notice also that the 
final central value of the NMC scalar field seems to de- 
pend logarithmically on the central density of the initial 
boson star <r(0) (in Fig. [3] we use a log scale for <r(0) in 
order to sec this logarithmic dependence). 



C. Gravitational collapse and black hole formation 

Since the scalarization phenomenon modifies the com- 
pactness of the star, one would expect that the threshold 
to the U-branch associated with the scalarized stars also 
changes depending on the value of £. For instance, the 
initial unscalarized configuration with er(0) ~ 0.266 and 
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FIG. 5: Scalar charge (top panel), and central absolute value 
of the NMC scalar field (bottom panel), for stable boson star 
configurations as a function of the initial central density of 
the boson field <r(0) and different values of £. Notice the 
logarithmic scale in er(0). 



£ = 1, which evolved into a quasi-stationary scalarized 
state as depicted in Fig. [2] will collapse into a black hole 
if we take instead £ > 10. This means that the criti- 
cal cr cr it(0) seems to decrease when compared to the GR 
case in a way that depends on the value of the param- 
eter £. In particular, we have found that for £ = 50, 
the last scalarized boson star configuration that is sta- 
ble against black hole formation corresponds to a central 
value er(0) ~ 0.195. This value of cr(0) has to be con- 
trasted with the value <r(0) ~ 0.272 for £ = (the GR 
case) , which is associated with the usual maximum mass 
configuration. In summary, depending on the value of 
£, there are two kinds of instabilities for boson stars in 
STT: On the one hand, there is an instability that takes 
the star into a stable scalarized state and, on the other 
hand, an instability that causes a scalarized star to col- 
lapse to a black hole. Presumably, an unscalarized boson 
star that collapses directly into a black hole may also 
reach a transient (possibly very brief) state of scalariza- 
tion before collapsing. 

We have studied the dynamical collapse of boson stars 
in STT toward a black hole. In order to be sure that a 
black hole is formed, we look for the appearance of an 
apparent horizon during the simulations. As expected, 
for unstable configurations an apparent horizon appears 
suddenly after some evolution. Its then grows for 
some time as more matter is accreted, until it finally 
settles once all the initial matter has either fallen into 
the black hole or has been radiated away. Figure [6] shows 
the apparent horizon mass (defined in terms of its area 
A as A/ah = yj A/16ir) as a function of time, for the first 
unstable boson star configuration for £ = {1, 50, 250, 500} 
and cr(0) = {0.268,0.195,0.231,0.249}. These values of 
the central density are above the threshold given by Table 
[U The asymptotic value of the NMC field </>o is dependent 
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FIG. 6: Apparent horizon mass of the first scalarized unstable 
boson star configurations for different values of the parameter 
£. When £ = 1 the configuration is essentially that of GR. For 
£ = {50, 250, 500} the scalarized boson star collapses toward 
a black hole of the same mass but at much earlier times. 



of £ through Eq. (|7.4p . Notice that as £ increases, the 
black hole formation time first decreases up to a certain 
value of £, after which the time of collapse increases again 
for larger £. 

Figure [7] shows some snapshots of the evolution of the 
NMC scalar field for a case where <r(0) = 0.23 and £ = 50, 
and an asymptotic value of <fio = —9.25 x 10~ n . A similar 
initial configuration with £ = (the GR case) is stable, 
while configurations with £ < 50 evolve into a stable 
scalarized state. 




FIG. 7: Snapshots of the evolution of the NMC scalar field 
for an initial configuration with rj(0) = 0.230 and £ = 50. 
The star eventually collapses into a black hole. The same 
initial configuration with <r(0) = 0.230, but taking £ = (the 
GR case), is stable. The vertical lines in the bottom panels 
indicate the location of the apparent horizon. 



D. Scalar radiation 

Gravitational radiation can carry energy away from an 
isolated system, and it also encodes important informa- 
tion about the physical properties of the system itself. 
In [59l | . Harada et. al have performed a numerical study 
of the scalar gravitational radiation emitted during an 
Oppenheimer-Snyder collapse in STTs in terms of the 
initial parameters, such as the initial radius and mass of 
the dust. 
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FIG. 8: Evolution of the Misner-Sharp mass for a stable bo- 
son star with central density cr(0) ~ 0.141. The variation of 
the Misner-Sharp mass with £ = 1 is insignificant while the 
variation of the mass with £ = 500 is around 20% (see table[T|. 
Notice however that the peak in the plots is not physical since 
uims acquires a physical meaning only in the limit r — > 00. 



In order to study the emission of scalar gravitational 
radiation in our boson star configurations, we will start 
by considering the reduction of the Misner-Sharp mass 
during the scalarization process. Figure [8] displays the 
evolution of the mass function rriMs(f) fo r an initial bo- 
son star configuration with central density c(0) ~ 0.141, 
both for £ = 1 and £ = 500. It is clear that the con- 
figuration with £ = 1 reaches rapidly a quasi-stationary 
state where the variation of the Misner-Sharp mass can 
be ignored (it is less than 1%). On the other hand, the 
evolution for the configuration with £ = 500 is quite dif- 
ferent. First, one can notice a distortion moving outward 
for which the mass function even becomes negative. The 
position of this distortion can be seen to coincide with an 
outward moving pulse of scalar field <f>. One should not 
worry about the fact that the mass function is negative 
since, as we have said before, one can only interpret this 
as a mass in the vacuum regions. More importantly, once 
this pulse has moved away, we can see a very clear reduc- 
tion in the mass function. At the end of the simulation, 
the final reduction in the mass is of around 20% (see ta- 
ble [TJ . We should also mention the fact that in this case 
the metric components and the scalar field <j> continue to 
oscillate at late times around a fixed configuration (these 
oscillations arc quite independent of the resolution of the 
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numerical evolution). 

One could naively think that this reduction in the mass 
is entirely produced by the emission of scalar gravita- 
tional radiation. However, at this point it is difficult to 
separate between the total amount of energy carried away 
by the scalar field <j> and that carried by the scalar grav- 
itational radiation itself. The latter of course arises due 
to the NMC between the scalar field and the curvature, 
but everything is mixed in the flux of energy as given by 
Eq. (|2.20[) . For instance, in the SS case with <fi = some 



scalar field is still radiated away during the transition to 
the quasi-stationary scalarized state (or during the col- 
lapse to a black hole) , even though in that case we do not 
expect scalar gravitational radiation since F '(^>)^ = o = 0. 
In order to have an unambiguous quantification of the 
amount of energy emitted in the form of scalar gravita- 
tional radiation one would need to go to second order 
perturbation theory and compute the equivalent of the 
Isaacson energy-momentum tensor in STT. We will leave 
this computation for a future work. 
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TABLE I: Initial and final values of the Misner-Sharp mass and integrated energy flux of the last stable boson star configuration 
in SST as measured by an Euler observer, for different values of £ which, according to Eq. (|7.4|l . correspond to different values 

of </)(). 



Figure [9] shows the final Misner-Sharp mass for a se- 
quence of boson star configurations both in GR (£ = 0) , 
and for STT with different values of £. The configura- 
tions presented here correspond to the S-branch of a sin- 
gle boson star for which the system has reached a station- 
ary state. The stationary unstable configurations (i.e. 
the U-branch) would continue on the right after the last 
plotted point of the curves. We have performed several 
numerical simulations using central amplitudes beyond 
the new S-branch in STT, and we have confirmed that 
the stars collapse into a black hole. 

Table U summarizes the initial and final states of the 
last stable boson star configurations found for each curve 
of Fig. IH1 Notice that for many of the cases studied here, 
at late times the system reaches a oscillating state, with 
oscillations that seem to decay very slowly, which makes 
it difficult to determine the final value of the Misner- 
Sharp mass. So, in order to provide an error estimation 
on the values reported on it, we consider a time average 
between the maximum and minimum value of the "final" 
mass. We also include the value of the integrated energy 
flux as measured by an Eulerian observer. Notice that 
the difference between the initial and final masses agrees 
very well with the energy flux for almost all the values 
of £, except for £ = 1. In this case the oscillations of the 
system do not allow us to get an accurate value for the 
final mass of the system. 

The evolution of the stars considered above proceeds 
in general as follows: After the initial burst of scalar radi- 
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FIG. 9: Misner-Sharp mass for a sequence of stable stationary 
scalarized boson stars for different values of £. For reference 
we also show the GR configurations (£ = 0). The last point on 
every curve corresponds to the last stable configuration found. 
To the right of that point all configurations are unstable. 



ation due to the scalarization process, the system settles 
down and reaches a state with very long-lived oscillations 
with a characteristic frequency. In the bottom panel of 
Figure [10] we show the main frequencies of this oscilla- 
tions for the stable configurations presented in Fig. [H] We 
obtain these frequencies by performing a discrete Fourier 
transform in time (at late times) of the non-minimal 
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scalar field at a fixed radius r = 120, and looking for the 
largest peak. One would expect that the scalar gravita- 
tional waves should have the same frequency [5!| . Notice 
that the main frequency of the system seems to be inde- 
pendent of the parameter £. The top panel of the Figure 
shows the actual Fourier transform (the power spectrum) 
for the case cj>o = 0.07 where one can see that there is in- 
deed a very clear peak. 
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FIG. 10: Top panel: Fourier transform in time of the NMC 
scalar field evaluated at r = 120 for different values of £ at late 
times (after the scalarization process). The initial configura- 
tion corresponds to a boson star with a(0) = 0.07. Bottom 
panel: Main frequencies (peaks in the power spectrum) for 
different values of cr(0). Notice that these frequencies seems 
to be independent of the parameter £. 



VIII. CONCLUSIONS 

Boson stars are stable self-gravitational configurations 
of a complex massive scalar field that evolves according 
to the Klein-Gordon equation. The expected mass of 
these type of stars typically varies between the mass of 
an asteroid and a few solar masses, depending on the 
mass parameter of the bosonic scalar field, some solar 
masses (38|- Though hypothetical, boson stars are sim- 
ple models that can be used to understand the corrections 
to general relativity proposed by alternative theories of 
gravity. Any alternative theory of gravity has to be tested 
against observations, both at the scale of the solar system 
and at the cosmological scale. Scalar-tensor theories of 
gravity, where a (real) scalar field is non-minimally cou- 
pled to gravity, are interesting generalizations of general 
relativity that have so far not been ruled out by observa- 
tions. 

In this paper we have used a scalar tensor theory of 
gravity to study the evolution of distinct families of single 
boson stars parameterized by the central density of the 
star and the parameter £ that controls the non-minimal 
coupling of the theory. We have focused on the transi- 
tion to a scalarized state using a fully relativistic code in 



spherical symmetry. We have found that, just at it hap- 
pens with neutron stars, boson stars can also undergo 
both spontaneous and induced scalarization, the latter 
case linked directly with the emission of scalar (monopo- 
lar) gravitational waves. 

Our numerical experiments show that the final mag- 
nitude of the non-minimally coupled scalar field seems 
to depend logarithmically on the central density of the 
boson star. We have also found that there is a critical 
value of the non-minimal coupling parameter £, corre- 
sponding to about to £ ~ 50, that maximizes the scalar- 
ization (i.e. the amount of final scalar charge). On the 
other hand, the maximum reduction of the initial mass 
of the boson star (about 20%), due to energy radiated 
to infinity during the scalarization process, was obtained 
for the maximum value of the parameter £ considered in 
our evolutions (£ = 500). We summarized our results in 
Table [IJ Each configuration considered there corresponds 
to the critical density that separates the stable and un- 
stable branches of a single boson star in scalar tensor 
theories for different values of £. Figure shows the sta- 
ble branches as a function of this parameter, once the 
system has reached final quasi-stationary state. It is ev- 
ident that, whereas a small value of £ leads to results 
that are very close to those of general relativity (£ = 0), 
larger values lead to important deviations that in princi- 
ple could be measured. 

At this point it is important to stress the fact that in all 
the evolutions presented here there are two parameters 
associated with the specific form of the scalar tensor the- 
ory: The free parameter £ used in the expression for the 
non-minimal coupling function F(<fi) = l + 87rGo£</> 2 , and 
the asymptotic (cosmological) value of the scalar field <fio 
that in our evolutions is chosen as the maximum value al- 
lowed by the constraints imposed by the Cassini probe on 
the effective Brans-Dicke parameter [given by Eq. (|7.3I) ]. 
This ensures that all our results satisfy the bounds im- 
posed by the Solar System experiments. 

Finally, we have taken a first step in trying to charac- 
terize the monopolar gravitational waves emitted during 
the scalarization process. One would expect that the 
magnitude of this monopolar radiation will be propor- 
tional to the product of the scalar charge and the square 
magnitude of the frequency of the radiation. We have 
found that this frequency seems to be independent of the 
coupling parameter £. 



Appendix A: Characteristic Variables for the BSSN 
Formulation in STT 



Using the slicing condition (|2.26p . and assuming that 
the shift vector is an a priori given function of the coor- 
dinates, Salgado et. al. have presented the characteristic 
decomposition of the BSSN formulation in the STT con- 
text [271 ] . They have shown that this formulation leads 
to a well-posed Cauchy problem in the Jordan frame. In 
this appendix we show that the algorithm used for deal- 
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ing with the regularization of the origin, in which one 
introduces a set of auxiliary variables to impose all the 
required regularity conditions [13, E3, [HI , does not spoil 
those results. 

Let us start by considering the spatial metric in spher- 
ical coordinates written in the form 



A(t, r)dr 2 + r 2 B(t, r)dVt 2 (Al) 
e 4x [a[t, r)dr 2 + r 2 b(t, r)dn 2 ] , (A2) 



where e x is the conformal factor and dil 2 the standard 
solid angle element. In order to construct a fully first 
order and regularized BSSN system, one first needs to 
introduce the following set of variables 



da = d r \na , 



d r X, 



d a = - d r \na . 
T = d r x , 



A r = - 
a 



dra_drb_ 2rX 

2a b 



(A3) 
(A4) 

(A5) 



where A is defined below in Eq. (|A15[) . It turns out that, 
with the above variables, the principal part of the BSSN 
formulation in the STT context, coupled with the slicing 
condition (|2.26|) . is given by 
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(A7) 
(A8) 
(A9) 
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(All) 
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where we have defined 8q — d t + (3 r d r . 

Notice that Eqs. (|A13IIA14|) are the spherical reduction 
of the Eqs. of motion (|2. 14112 .15)) . Moreover, the quanti- 
ties A and A\ arc auxiliary variables needed in order to 
impose regularity conditions at the origin, and which are 
defined by 



A 



1 



(1 - ajb) 



A, 



1 



(A a - A b ) . (A15) 



Also, the equation of motion for A r has been modified 
by adding a multiple of the momentum constraint which 
is controlled by the r] parameter. In the following, we 
consider, for simplicity, that r\ = 2 (see (20j). 

Taking into account the fact that an evolution system 
of the form 



dtu = v\d r v , dtv = v^drU , 



(A16) 



has the characteristic variables w± = it=F \/ (^2/^1) v with 
characteristic speeds ^2, it is easy to show that the 

subsystem (|A13IIA14)) has the characteristic decomposi- 
tion 



Uf =Q r ±e 2x v^n, 



with speeds of propagation 



A4 = P r T > 



-2x 



(A17) 



(A18) 



On the other hand, the characteristic variables related 
with the choice of slicing, and associated with the sub- 
system (|A6[) and (|A9|) . turn out to be 



U: 



± 



d a ± ae 2x yj af B M K 
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x {Q,.(/bm - 0) - e 2 Va/BM(l - 0)U}, (A19) 
with characteristic speeds 
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Furthermore, the following combinations 



Ui = d x ± 



f {Q r (l - 2/bm + 0)- VEe 2x {l - 6)U} 
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a ( 1 „ r 2 d 
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, (A21) 



provide two more characteristic variables with associated 
speeds 

M± = fi r T a e -^- . (A22) 
v a 

Finally, the following variables 

a J V o/bm/ oa/BM 
V = _ aA r + f'Qr{2f'f BM -0) + d^ + 2T ^ (A24) 
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are the remaining characteristic variables with speed 



(A25) 
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Appendix B: Constraint preserving boundary 
conditions 



To reduce the influence of the numerical boundary 
on the dynamics of the system we have considered con- 
straint preserving boundary conditions (CPBC), in the 
sense that small violations of the constraints introduced 
by spurious reflections from the boundary converge away 
with the resolution. This is not only motivated by the 
requirement of having a well-posed initial boundary prob- 
lem, but mainly because for long evolutions such as those 
considered here, constraint violating modes may propa- 
gate inside the domain contaminating the interior solu- 
tion, which could lead to an incorrect physical interpre- 
tation |62l ]. 

In order to impose our CPBCs we introduce informa- 
tion from the characteristic variables at the boundary 
through the following algorithm. 

1. We construct numerically the characteristic outgo- 
ing modes U+ using the dynamical variables, such 
as gauge, metric components and </>, at the bound- 
ary. For instance, one can reconstruct the outgoing 
mode U+ at the boundary by using one-sided dif- 
ference scheme for II and Q r . 

2. In order to reconstruct the incoming scalar and 
gauge fields, Uf_ and [/^ ause , we assume that <j> and 
a behave at the boundary as outgoing spherical 
waves of the form: 



u(t, r) =u + - f(r - t) . 
r 



(Bl) 



with uo their corresponding asymptotic values. By 
taking time and space derivatives one can then re- 
construct the corresponding incoming characteris- 
tic fields. For example, for U_ we find (assuming 
that <fi = asymptotically): 



,-2x 



Ul = — 



' ar 



-<t>- 



(B2) 



It is interesting to notice that the incoming fields 
arc not zero, as one could naively expect (just set- 
ting them to zero introduces quite large reflections). 

3. This leaves us with the incoming field U_. This 
is where one can impose the constraints, as one 
can show that the spatial derivative d r U* can be 
written as a combination of the constraints plus a 
term that contains no derivatives of the dynamical 
variables. Asking then for the constraints to vanish 
at the boundary allows us to evaluate d r U* directly 
at the boundary. And once we have d r U* at the 
boundary we can use it, together with the values 
of at the nearby points, to solve for U_ at the 
boundary by simple finite differencing. 



4. Finally, we recover the dynamical variables using 
both the incoming and outgoing modes. Notice 
that it is not necessary to recover all the variables, 
since many of them are in fact not independent. 



The above algorithm allows us to impose CPBCs for 
the spherical reduction of the BSSN formulation. As an 
example, in Figure [TT] we show results from an evolution 
of a single boson star using forth order spatial differ- 
ences and a forth order Runge-Kutta for the integration 
in time. The outer boundary is locates at r out = 240. 
The figure shows the Li norm of the violation of the 
Hamiltonian and momentum constraints for a boson star 
with cr(0) =0.1, and for two different resolutions. No- 
tice that while the light crossing time of the numerical 
grid is t ~ 240, the scheme remains fourth order conver- 
gent for much longer times. The boundary conditions do 
a spurious reflection whose magnitude is less than 10~ 9 . 
This effect is only evident in the momentum constraint at 
t = 240 and t = 480. The magnitude of the violation of 
the Hamiltonian constraint in the interior is much larger 
than those reflections. 
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FIG. 11: Violation of the constraints for a boson star with 
cr(0) = 0.1. The Li norm of the Hamiltonian (top panel) and 
the momentum constraints (bottom panel) is plotted for two 
different resolutions as a function of the time. The effect of the 
outer boundary conditions is visible only in the momentum 
constraint at times t = 240 and t = 480. 
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i,j, k, ■ ■ ■ will denote 3-dimensional spatial quantities. 



